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I. Introduction 



The sine kernel 



sm i - y ) 

(N 1 K{x,y):=— 

7t{x - y) 

arises in many areas of mathematics and mathematical physics. There is an extensive 
literature on the asymptotics of the eigenvalues of K s , the operator with this kernel 
on an interval of length s, as s — > oo, for example f|, H, H, @, and asymptotic fomulas 
of various kinds were obtained. Some of these derivations were rigorous, others were 
more heuristic. 

; The Fredholm determinant of the kernel is of particular interest. In the bulk 

scaling limit of the Gaussian Unitary Ensemble of Hermitian matrices, the probability 
that an interval of length s contains no eigenvalues is equal to 



det (/ - K s 



The asymptotics of Fredholm determinants of convolution kernels k(x — y) have a 
^ ■ long history. (The history of their discrete analogue, Toeplitz determinants, is even 



longer, beginning with the 1915 paper [O of G. Szego.) If the Fourier transform k 



of k is smooth and less than 1, and if k satisfies some other conditions, then for the 
corresponding operator K s one has as s — > oo 

log det (J - K 8 ) = ci s + c 2 + o(l), 



where c\ and c 2 are explicitly determined constants [10|. If k < 1 and is smooth 
except for jump discontinuities then the result becomes 0] 

log det (I - Kg) = ci s + c 2 log s + c 3 + o(l). 

Even for some cases where k < 1 is violated at finitely many points a similar relation 
holds 0. But for the sine kernel k is the characteristic function of the interval [—1, 1] 
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and so the condition is violated on a set of positive measure and the situation is 
considerably more complicated. 

The asymptotics of the determinant for the sine kernel were probably first inves- 
tigated by Dyson ||. He derived the asymptotic formula 

det (/ - K s ) ~ 2 1 / 3 e 3 ^- 1 ) s- 1 ^ e~ s2 '* (1) 

by applying a scaling argument to a known asymptotic formula for Toeplitz determi- 
nants fl6| , and then used inverse scattering techniques to complete this to an asymp- 



totic expansion. In 1980 Jimbo et al. || showed that the logarithmic derivative of 
this determinant satisfies a second-order differential equation which is reducible to a 
Painleve V equation. The asymptotic expansion can also be obtained, starting from 
([!]), by substituting a formal expansion into the equation and successively solving for 
the coefficients. 

Another quantity of interest from the random matrix point of view is 



tr K s (I - K s 



which in the same model equals the probability that an interval of length s contains 
precisely one eigenvalue, divided by the probability that it contains none. (There are 
similar but more complicated formulas for any finite number of eigenvalues jnj].) The 
asymptotics of this were derived in 0] where 

trK.V-K.)- 1 *-*, (2) 
Zy 2irs 

was obtained by scaling an anologous result for Toeplitz matrices, which in turn was 
proved by exploiting its connection with orthogonal polynomials, whose asymptotics 
were known. This formula was then extended to a complete asymptotic expansion by 
using the differential equation. 

None of these asymptotic results were rigorously proved by the methods described. 
In [p|j , by studying the asymptotics of a continuous analogue of orthogonal polyno- 
mials, we were able to give a proof of (||D and of the first-order asymptotics in (|l|), 

log det (I - K s ) s 2 /8, 

actually the slightly stronger 

log det (J - K s ) = -s/A + 0(1). (3) 

as 

As far as we know nothing beyond this has actually been proved up to now, and 
refinements of these are not the subject of this paper. The subject is the extension 
of this work from the one interval case to several intervals. Thus a single interval 
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of length s is replaced by sJ, where J is a fixed union of intervals. Even these 
asymptotics are quite elaborate and involve hyperelliptic integrals and the Jacobi 
inversion problem for Abelian integrals. 

We continue to denote our operator by K s , and it is understood that it now acts 
on sJ. We shall show that as s — > oo 

- — log det (J - K 8 ) = ci s + 02(a) + o(l), (4) 
as 

where c\ is a constant and 02(5) is a certain bounded oscillatory function of s, and 
that 

tr K s (I - K s y l = Oc 4 (s) + GOT*)), (5) 

where C3 is another constant and c^s) is another bounded oscillatory function. The 
constants C\ and C3 are explicitly computable but the determination of the functions 
C2(s) and 04(5) requires the solution of a Jacobi inversion problem. (At least this is the 
case with this method. It is possible that a different approach might lead to simpler 
representations.) If J consistes of m intervals then there is a curve, parametrized by 
s, in the m — 1-torus T m , and two real- valued functions defined on T m . These 
functions, when restricted to the curve, are Oi(s) and c^(s). When m = 2 they are 
periodic, and we can write down an integral representation for the period. 

Here is how we obtain the asymptotics. For the sequence of monic orthogonal 
polynomials P n (z) associated with a weight function w on the unit circle, the square 
of the L 2 (w)-norm of P n is the ratio of two consecutive Toeplitz determinants asso- 
ciated with w. The asymptotics of P n gives information on the asymptotics of this 
ratio of determinants. Analogously the asymptotics of the continuous analogue of 
these polynomials, which we shall denote by R(x), determine the asymptotics of the 
logarithmic derivative of the Fredholm determinant. The underlying weight func- 
tion for R(x) will be the characteristic function of the complement of J in R, so we 
have the continuous analogue of the polynomials orthogonal on a union of arcs. The 
asymptotics of orthogonal polynomials such as these were obtained in |T5[ and the 



ideas of this paper were used in [|17]] and are used here. We find an entire function 



which comes close to satisfying the characteristic property of R given in (y) below by 
first finding another function (which is not entire, or even single-valued) which has 
this property exactly (the function h s in Lemma 4), and then we approximate this 
by an entire function. The entire function we find in this way will necessarily be a 
good approximation to R. This gives a weaker form of @ at first. It is strengthened 
by exploiting its relation to a certain extremal problem. Some details will be omitted 
here since the complete details for the analogous orthogonal polynomial case can be 
found in |T5|| . 

We obtain (|5|) by using a representation, which is essentially contained in |9[], for 
the resolvent kernel of K s in terms of R(x). There are analogous representations for 
a large class of kernels like the sine kernel [0, H . 
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In || it was shown that for general J the Fredholm determinant is governed by a 
system of partial differential equations with the end-points of J as the independent 
variables. The dependent variables are the values of the function R(x) at the end- 
points of J, and the logarithmic derivative of the Fredholm determinant has a simple 
representation in terms of these values. We do not use this representation but rather 
the one alluded to above (formula (|7|) below), which makes its asymptotics, at least 
by this method, more accessible. Nevertheless, our results may give a hint of what 
the asymptotics of the solutions of the system of equations might involve. 

The author wishes to thank Craig Tracy for introducing him to the subject of 
random matrices, and in particular to the problem of asymptotics for the several 
interval case. 

II. The Resolvent Kernel and the Function R(x) 

We shall denote the Fourier transform by a circumflex, as usual, 

/■CO 

/(*):= / e ixt f{t)dt, 

J — oo 

and write £ s for the space of Fourier transforms of functions in La(— s, s). This 
consists of entire functions which are 0(e s ' Imx ') and whose restriction to the real line 
R belongs to L 2 . We also write E for the complement of J in R and Q for the 
complement of E in C. 

Lemma 1. There exists a unique function R(x) G e~ %sx + S s such that 

/ JWx) - e~ isx ) g(x) dx= [ g(x) e isx dx for all g E £ s . (6) 

JE J J 

In terms of this function we have the representations 

-4-logdet(I-iT s ) = -|J|+- / \R(x) -e~ isx \ 2 dx } (7) 

ds IX IX JE 

where | J\ is the measure of J, and 



tr K s (I - Ks)- 1 = Im / R\x) R(x) dx. (8) 



Remark. The set e~ lsx + £ s is analogous to the set of monic polynomials of a given 
degree, and @ implies that R(x) is orthogonal on E to the Fourier transforms of a 
dense set of functions in L 2 {— s, s), those which are smooth and vanish at — s. (The 
integral expressing the orthogonality is a principal value integral at infinity.) Thus 
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R(x) is analogous to a monic polynomial of a family orthogonal on a union of circular 
arcs. The identity fl7|) is the analogue of the fact that the square of the norm of the 
monic polynomial is the ratio of two consecutive determinants. 

Proof. Write 

sin s(x-y) 

K s {x,y) := r— , 

tt(x - y) 

and denote by K s the operator on all of R with kernel K s (x, y) Xj(y)- Our determi- 
nant and trace are the same for this kernel as for the other. We express the resolvent 
kernel R s (x,y) for this operator in terms of the functions 

R±{x) :={I-K s y l e ±isx 

introduced in M. If M denotes multiplication by x then the commutator [K 8 , M] 
has kernel 

\ (Jsx -isy _ „-isx isy\ 

2irV '' 
It follows from this upon left- and right- multiplying by (I — K & )~ 1 that for x,y G J 

R+jx) R_(y)-R_(x)R + (y) 

R s (x J y) = r , 

2ni(x — y) 

and in particular 

n , \ R'Jx) R-(x) - R'Jx) R+(x) 
R s {x,x) = — 



2ni 

We write R(x) for R_(x) and observe that, since the kernel of K s is real, R+(x) 



R(x). Thus 



R(x) R(y) - R(x) R(y) 
Ks ^ y >- 2m{x-y) 



R s (x, x) = Im R'(x) R(x). 

This establishes (H). Next, since the kernel of K' a (differentiation is with respect to 
s) is 

J_( e *-(*-v) +e -*.(*-i/)) Xj ( y ) j 
the kernel of (J — K^ 1 K' s equals 

^(R+(x) e- is y + R.(x) e ls y)Xj(y) = -Re R{x) e^Xj(y). 

Z7T 7T 



Thus, since 



4- log det (/ - K s ) = -tr (/ - K s )- 1 K 1 , 
as 



we obtain 

- jj- log det (I - K s ) = i Re jf fl(ar) e isa: dx. (9) 

If we apply the identity @, which will be proved momentarily, to the function g{x) = 
R(x) — e~ tsx it gives 

/ \R(x) - e~ isx \ 2 dx = [ R(x) e isx dx - \J\, 

JE J J 

and (0) follows from @. 

To obtain the first assertion of the lemma observe first that the difference between 
any two functions R satisfying belongs to £ s and is orthogonal on E to £ s , and 
so it vanishes on E. Since £ s consists of entire functions, this difference must be 0. 
To show that our fuction R(x) = satisfies we write its definition in the 

equivalent form 

R(xj- J j K s (x,y)R&dy = e iSX 

and observe that K s (x,y) is the kernel of the projection operator from /^(R) to £ s . 
Hence g G £ s implies 

/ g(x) K s (x,y)dx = g(y). 

JR. 

Thus if also g G Li(R) then multiplying both sides of the previous identity by g(x) and 
integrating over R give (|6]) in this case. The extension to general g is straightforward. 
□ 



II. Green Functions, Neumann functions and 
a Reproducing Function for E 



In [171, where J 



was to show first that 
values on E of the function 



— 1,1], E = (oo, — 1] U [l,oo) and Q = C\E, the procedure 
H) holds if R(x) is replaced by the sum of the two limiting 



h(x) 



s\/ x 2 — l 



X + a/x 2 — 1 

\ 2 V / ^T" 



(10) 



defined in Q with appropriate branches of the square roots. Notice that R(x) and 
h(x) have the same asymptotic behavior as x — >• oo in Q. One approximates h(x) 
by an entire function, for which (^j holds approximately, and then deduces that this 
function must be a good approximation to R(x). 
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The function \/x 2 — 1 in (|T0D was chosen to be that branch of the square root 
which is analytic in Q and asmptotically equal to x as x — > oo in the upper half- 
plane. Its imaginary part vanishes on E and equals |Im x\ + 0(1) in Q, and so is a 
kind of Green function for Q. We begin by constructing its analogue for general 

m 

J= U K A]- 

fc=i 

Let 

fc=l 



and let yq(x) denote that branch which is analytic in f2 and asymptotic to x m as 
x — > oo in the upper half of Q. This has purely real limiting values on E. Let p(x) 
be the monic m'th degree polynomial determined by the m equations 

/, *_P^L dy = 0) (* = l,...,m). (11) 
Then define 

G (x):= [ X ^Ldy, (xefl). (12) 

Jqi yg(y) 

This is multiple- valued when m > 1 because of the presence of the intervals aik+i] 
(A; < m — 1) one might integrate around and get a nonzero result. If Ck is a curve 
going around this interval and none of the others then the periods 

A Ck G =i -^Ldy (13) 

J lPk,a k+1 ] Jq(y) 

are purely real, so that Im G(x) is a single-valued harmonic function in Q. The 
conditions QTTJ) guarantee that Im G(x) — > clS X ^ 
For some real constants a[ and 02 we have 

P £L = l + ^ + ^ + 0(l) asxtoo. 
q[x) xx x 



Here the notation x f oo means that a; — > oo in the upper half-plane of Q. (Similarly 
x I oo will mean that x —>■ oo in the lower half-plane of Q.) The integral 



oo+Oi 
-oo+Oi 



P{X) -l\dx 



'q(x) 
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equals —iria^. But it follows from (|IID that it is also real. Hence a[ = and so 

P (X)^ = 1 + 02 + 0( J_ } asxT(X)> (14) 



It follows readily that 



G fa;) = x + 




X X A 



where a±, the expression in brackets, is purely real. We now define 

G(x) := Go(x) — a±. 

G(x) has the same Green function characteristics as Go(x) and has the behavior 

G(x) = x + 0{— ) as x | oo. (15) 



Since Go is purely imaginary on J, we have Go fx) = — Go(x), and so from the above 
we deduce 

G{x) = —x — 2 di H — - + 0(— ) as x j oo. (16) 

The function e lsG ^ is multiple-valued and analytic in f2 and has single-valued 
absolute value. Such a function can be thought of as a section of a holomorphic line 
bundle over Q, but instead we proceed as follows. Denote by 7i the set of all (single- 
valued) analytic functions / in the complement of (— oo, cti] U [j3i, oo) in C for which 
the limits 

u + xe [a k ,(3 k ], (k = 2 r ..,m) (17) 

e^0+ J[X— IE) 

are constants of absolute value 1. Each such function continues to a multiple- valued 
function in Q with single-valued absolute value. Two functions in 7i are said to 
belong to the same class if the corresponding limits ( |TTD for them are the same. 
A class is denoted by V and the set of classes can be identified with the m — 1- 
torus. We denote by Ti-(T) the functions of class T. Occasionally we shall allow 
our functions to have poles. An example is the exponential Green function $(x, xo) 
described below. We shall write our classes additively, and so each will be an element 
of T m_1 := (R/2vrZ) m - 1 . If T is the class of e iG(x ^ then S r is the class of e isG(x \ 
This is our curve in T m_1 . 

Next we have to find an analogue of the second factor in ([TDD . There will be 
many, one for each class sr , and so we will find one for a general class T. We begin 
by recalling the characteristic properties of the Green and Neumann functions for 
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Q. Given xq, Green's function with pole at Xq, denoted by g(x, Xq), is harmonic for 
x G Q, except at Xq where it has the singularity log \x — £o| _1 ? and has limit as 
x — > E. Denoting by g(x, xo) its (multiple-valued) harmonic conjugate, we define the 
exponential Green function by 

$(x,x ) := e 9{x ' Xo)+i ' 9{x ' Xo) . 

Neumann's function N(x,xq,xi) is harmonic except at x = xq and x = x±, where it 
has the singularities log | (x — Xi)/(x — xo) | , and on E 

dN(x,x ,xi) _ 
dn x 

Here d/dn x denotes the two derivatives normal to E. We set 

^{x^o^t) := e m^o,x 1 )+iN(x,x ,x 1 )_ 

This has a simple pole at x = Xq and a simple zero at x = x\, unless the two are 
equal, in which case the function equals 1. The function $(x, x ) is multiple-valued 
with single-valued absolute value, and is determined only up to an arbitrary constant 
factor of absolute value 1. The function fy(x, xo, xi) is single-valued but determined 
only up to an arbitrary nonzero constant factor. The analogue of the second factor 
in ( |10D will be built out of these functions. 

Lemma 2. The function 1 + p(x)/q(x)^ is nonzero (i.e., has nonzero limits) on E 
and has m — 1 zeros in Q. 

Proof. We know that 

Im (x + G{x)) = on E. (18) 
From (|15D and ([16|) we deduce that 

lim^oolm (x + G(x) > 0. 

Therefore 

Im (x + G(x)) > in fi. (19) 

If 

l + ^=L = (x + G(x))' 
^Jq(x) 

had limit zero at some point of E this would contradict QT8| ) and ( |T9| ) since at a 
critical point an analytic function maps a local half-disc to a full disc. 

For the second statement let Ck be the contours in ( |i3|) and C = C- U C + , where 
C- starts at — oo — Oi, goes around oti, and ends at — oo + Oi, while C + starts at 
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oo + Oi, goes around f3 m , and ends at oo — Oi. All the contours are described so that 
the parts below E are traversed to the right, the parts above E to the left. We find 
that 



A c arj 



p(x) 



0, A Cfc arg 1 



p(x) 



/q(x)J \ ^q(x) 

This and the argument principle establish the lemma. 



-2tt, (k = 1, 



m 



□ 



Denote the zeros of l+p(x)/q{x)z by x\ 



in 



, x m _i*. It follows from a discussion 
15] that for each class Y there are unique points xi, • • ■ , x m _i, and unique numbers 



, £ m -i each equal to ±1, such that 



m— 1 

IJ ^(x,x fc *,x fc ) > on £ 

fc=0 

(with an appropriate normalization of the \l/'s) and such that 

p(x 



(20) 



kAx) 



1 + 




m— 1 

[] ^(x,x k *,x k )$(x,x k ) 

k=0 



e r. 



(21) 



The x k are given by the solution of a certain Jacobi inversion problem. The details, 



which we shall not present, can be found in §6 of [15 |. 

One possibility, which occurs for certain T, is that all x k = x k * so the \l/'s don't 
appear at all. In general, though, they do appear and the zeros of 1 + p(x)/q(x)2 
are cancelled by the poles of the product of the \l/'s. The zero x k of ^(x,x k *,x k ) 
is cancelled by the pole of $(x,x k ) £k if e k = 1 and reinforced to a double zero if 
e k = —1. So there is no extra multiple- valuedness introduced by taking the square 
root. Observe, though, that we have not quite defined kr because of the nonuniqueness 
of the exponential Green and Neumann functions. We make them unambiguous by 
the requirements 



N(x, x k *, x k ), N(x, x k *, x k ) 



and then choose the square root so that kr(x 
Our eventual replacenent for h(x) will be 



g(x,Xk*,Xk) — * as x | oo, 
• 1 as x | oo. 



h s (x) := e - isG W k sro (x). 



(22) 



(23) 



Lemma 3. Suppose if> e 7i(r) has the behavior 



c 1 
ip(x) = — h O(-) as x 1 oo, 



1 



x 



x^ 



i/j(x) = O(-) as x I oo. 

X 
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Then 

kr(x) ip(x) \dx\ = —tt i c. 



E 



Remark 1. The conjugate- analytic function kp belongs to the class — T so that 
kp(x)ip(x) is single-valued in Q. The notation § E ---\dx\ indicates that E is not 
oriented but is taken twice, using the two limiting values of the function. 

Remark 2. It is because of this caracteristic property of kp, reminiscent of that of 
reproducing kernels, that we call it the "reproducing function" associated with T and 
Vl: integrating ip against hp yields the constant determining ^'s behavior as x f oo. 

Proof. Let jp(x) be the function obtained by replacing each by — Ek in (|21|). Then 

\jp(x)\ = \k r (x)\ onE (24) 

and, by 



pix^j 

jp(x) kp(x) is real- valued on E with the same sign as 1 H — , (25) 

lq{x) 



From the fact that Im (x + G(x)) is zero on E and positive in Q it follows that 

l + ^L = (x + G(x))' 



q(x) 

has positive limit on E from above and negative limit from below. Alternatively, if 
we denote the "upper" (resp. "lower") part of E by E + (resp. E~), then 

p{x) \ f 1 on E + 
sgn|l + ^l=^ 4 QnE -_ 



1 on E + 
-1 onE~, 




Hence, from (pl|) and (|25|) 

kpkp = kpjp x 



since both sides have the same absolute value and are positive on E. Dividing by kp 
we see that the integral in the statement of the lemma equals 

jp(x) i/j(x) dx, 

E 
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where now E is oriented so that E + is traversed to the right, E~ to the left. Observe 
that jr € 7i(—T) and so jr(x)ip(x) extends to a single- valued analytic function in 
Q. Moreover 

{0(x~ 2 ) as x I oo, 

(26) 
ci" 1 + 0(x~ 2 ) as x 1 oo 

(the first because 1 + p(x) / \J q(x) = 0(l/x) as x | oo). The integral is equal to the 
limit of the sum of the integrals over two semicircles, one in the upper half-plane and 
one in the lower, as their radii tend to infinity. By (^6]) the integral over the upper 
semi-circle tends to while the integral over the lower tends to — 7r ic. □ 

Now recall that we define h s (x) by (p3|) . 

Lemma 4. Suppose g(x) is single-valued and analytic in Q and satisfies 

e isx g(x) = - + Of— ) as x f oo, e~ isx g(x) = O(-) as x | oo. 

Then 

<j> h s (x) g(x) \dx\ = —nic. 



Proof. Write the intgral as 

<fk^e isG ^g(x)\dx\. 

J E 

The second factor belongs to class sr and by (|I3|) and fllED has the properties of 
ip(x) in the statement of Lemma 3. The assertion follows. □ 



III. Asymptotics of R(x) 



Define, for x 6 E, 



h±(x) := lim h s (x + is) 

e^0± 



We shall deduce from Lemma 4 and the first statement of Lemma 1 that R is well- 
approximated on E by h + + To do this we first replace h + (x) + /i_(x) by an 
entire function with the same general behavior as R and which is close to h + + h_ on 
E. Such a function will be e~ lsx + q(x), where 

— isx pia+oo p ls y h(n\ 

q(x):= e — e -^dy (27) 

2irt Jia-oo y — x 

with a < min (0, Im x). (We have dropped the subscript s in the notation for the 
function h s .) To prove that we do get a good approximation we derive an integral 
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representation for q(x) which requires deforming the path of integration in a rather 
involved way, as well as the analytic continuation of h(y). 

Recall that the exponential Green functions in ([H]) have absolute value 1 on £ 
and by (|20|) the the product of exponential Neumann functions is positive there. Thus 
both of these continue to the Riemann surface obtained by joining Q to a copy Q in 
the usual way, by identifying E ± in Q with E T in Q. Similarly G(y) continues into Q 
and its continuation has negative imaginary part there. The resulting continuation of 
h{y) is not single- valued because of the presence of branch points at the ct^ and (5^. 
We shall deform our path of integration to a system of contours, all lying in Q where 
Im G(y) > 0, and all lying in the upper half-plane where the factor e tsy is small. 

Before stating the main lemma we describe the contours Ai, (i = 1, ■ • • , 2m). If 
% — 2k — 1, then Ai begins at + ioo and goes down toward a^, it makes 1 1/2 
infinitesimal counterclockwise circuits of and then goes back up to + ioo. If 
i = 2k then Ai begins and ends at j3k + ioo but goes around (3k, when it gets to it, 
clockwise instead of counterclockwise. Except for the infinitesimal loops around the 
points of d J, these contours lie entirely in the upper half-plane of Q, and so have the 
sought properties. 

The function h(y) in the integrand below will be a continuation of the original 
h(y). (This will be clear during the course of the derivation.) It can be described as 
that continuation of h which, at the initial points of the contours A iy is the function 
obtained by continuing our original h from below E across it into the upper half- 
plane; on the rest of the Ai we take the continuation along Ai of this one. Observe 
that when our contours first cross the real line through E, as they do, h(y) returns 
to its original value since we will have undone the first continuation. 

Lemma 5. For x G E we have 

e~ tsx + q{x) = h+ (x) + h- fx) + e~ lsx V ^ '— \ KiJ> dy. (28) 

2m J a, y - x 



Proof. The path of integration in (p7[) may be deformed to the union of the contours 
C and Ck which occur in the proof of Lemma 2 (all described so that E~ is traversed 
to the right and E + to the left) plus an infinitely large semi-circle in the upper half- 
plane traversed clockwise. Since e lsy h(y) = 1 + 0(|y| _1 ) as y | oo the integral over 
the semi-circle equals — ni, and so 

e^q{x) = -i + -L/ < —^y±dy. (29) 

2 2m JcuCiU-c m -i y — x 

Suppose now, for definiteness, that x G (— oo, ax), and consider C_, the left part of 
C, which goes from — oo — Oz, around oi\ counterclockwise, and then to — oo + Oz. We 
move the upper part of C_ through (— oo, a%) into the lower half-plane and the lower 
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part of C- through (— oo, a.\) into the upper half-plane. The result is a contour which 
runs from cti + ioo to cti, around «i counterclockwise, and then down to a.i — ioo, 
with h(y) in the integrand replaced by its continuation. If we keep in mind how G(y) 
and kr continue we see that the continuation of h(y) from the lower half-plane to the 
upper satisfies 

e^%) = 0(^), 
while its continuation from the upper to the lower satisfies 

e isy h(y) = l + 0( T -r). 

\y\ 

This shows that the deformation is valid if we add to the integral the integral over an 
infinitely large quarter-circle in the lower half-plane, which is —iri/2, as well as the 
contributions from the pole at y — x, which is passed twice during the deformation. 
It follows from this discussion that if we replace C_ by this new contour, which we 
call Bi, then we must add 

~ + e isx h+(x)+e isx h_(x) 

to the right side of (p9|) . Next, the contour C\, which runs counterclockwise around 
[Pi, a 2 ], we deform by pushing its upper part down through (Pi, a 2 ) into the lower 
half-plane and its lower part up through (Pi, a 2 ) into the upper half-plane. Both 
continuations of e tsy h(y) are bounded at 00 and the result is that we can replace C\ 
by — B 2 + B 3 , where B 2 (resp. .83) starts at (3\ -Moo (resp. a 2 -Moo), loops around (3\ 
clockwise (resp. a 2 counterclockwise), and ends at (3\ — ioo (resp. a 2 — ioo). During 
this deformation no extra terms are picked up. We continue analogously with the 
remaining contours and then end with C + , the right part of C, deforming it as we 
did C_. [The reader is advised to draw a picture.] Another integral over an infinitely 
large quarter-circle is picked up during this last deformation. The result at this stage 
is 

e tsx q(x) = -l + e tsx h + (x) + e lsx h„(x) + V ^ '- — / —dy. (30) 

^ 2ni Jb, y — x 

Finally the lower parts of B\ and — B 2 can be joined into a single contour, and the part 
of this new contour which runs from a,\ to (3\ in the lower half-plane can be pushed 
up through («!, Pi) into the upper half-plane, resulting in the contour A\ — A 2 as in 
the statement of the lemma. This deformation is valid since G(x) continues to have 
negative imaginary part during it. The remaining contours are deformed similarly, 
and this completes the demonstration of the lemma. □ 
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Lemma 6. We have, as s — » oo, 

+ q( x ) — h + (x) — h_(x)\ 2 dx = 0(s~ 1 ). 



Proof. We must show that the norm in L2{E) of the sum in (p8|) is 0(s~5). Consider 
the integral over Ai, which can be written as a single integral paramaterized as 
y = cki + it, t > 0. (The integral over the infinitesimal loops around «i vanish.) For 
each 5 > the integral over i > 5 is bounded by a quantity which is exponentially 
small in s, uniformly in x, times 0((1 + M) -1 ), so the L 2 (E) norm of this part of 
the integral is exponentially small. In the neighborhood of t = 0, Im G(a\ + it) is 
asymptotically a negative constant times £2 (see the remark below), while 

\ksr (cxi+it\ =0(r 1 / 4 ). (31) 

(This last comes from the behavior of 1 +p/^/q-) Since on A\ 

1 



O 



\x — y\ \ \x — ai\ + 1 1 
we see that this part of the integral is bounded by a constant times 



t 1 / 4 (\x - Oii\ +t) 



dt (32) 



for some r] > and it is a simple exercise to show that the /^(R) norm of this is 
0(s~2). The other A{ are treated in a similar manner. □ 

Remark. Since the function k s r in (|31| ) depends on the parameter s we have to 
know that the estimate holds uniformly in s. The reason it does is that the Green 
and Neumann functions appearing in (^) are smooth, uniformly in their parameters, 
away from their poles in the Riemann surface ft U fl. (See, for example, Lemma 4.1 



of (T^j.) The local coordinate for the surface near x = ct\ is {x — ct\)^ and it follows 
that k s r and G(x) are smooth functions of (x — 0:1)2. This shows that ([31]) holds 
uniformly in s and also accounts for the yt in the exponent of the integrand in (B3). 



Lemma 7. We have, as s — ► 00 



/ \R(x) - h + (x) - h^(x)\ 2 dx = 0(s- 
Je 



Proof. The integrand equals 

(R(x) - h+(x) - h-(x)) (e- isx + q(x) - h+(x) - h_(x)) 
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+(R(x) - h+(x) - h-(x)) (R(x) - e~ isx - q(x)). 

Let us look at the last term. It follows from Lemmas 1 and 4 that the integral over 
E of the first factor times the complex conjugate of the Fourier transform of any 
smooth function in L/2(E) vanishes, and by continuity that Fourier transform does 
not have to be of a smooth function. Now R(x) — e~ lsx is the Fourier transform of a 
function in L,2(E), and so is q(x): it is easily seen from its definition ( p7[ ) that q(x) is 



an entire function of exponential type s and it follows from Lemma 6 that it belongs 
to L 2 (R). Thus the integral over E of the second term above vanishes. By Lemma 
6 and Schwarz's inequality the integral of the first term is 0(s~z) times the square 
root of the integral we started with. This establishes the lemma. □ 

This is all we shall need for @ and the reader interested only in these can go on 
to the next section. A few more estimates will be needed for the proof of (|5|). 

Lemma 8. (a) e~ lsx + q(x) = h(x) + 0(s~ 3 ^ 2 ) uniformly on compact subsets of fi. 
(b) For every e > there exists a neighborhood of d J in which q(x) = 0(e £S ). 

Proof. There is a modification of (Pq), proved in the same way, which holds for 
i6fi: the sum h + (x) + h-(x) is replaced by h(x) and the contours A4 are adjusted 
so as to avoid x (or the compact subset of Q in which they may lie), but they still lie 
in the upper half-plane of Q. This is easily seen to give assertion (a). It follows as a 
special case of this that if S is sufficiently small then q(x) = 0(e £S ) on the circular arc 

x = a k + Se l , — — < 9 < — . 



For the rest of the circle we use the continuation of Q2q), as originally stated, from 
x G E onto that arc and obtain a similar estimate. Since q is analytic inside the circle 
it is 0(e £S ) in the full disc. A similar argument applies, of course, to neighborhoods 
of the p k . □ 

Lemma 9. (a) R(x) = h(x) + 0(s~5 e s!mG ^) uniformly on compact subsets of the 
interior of J. 

(b) For every e > there exists a neighborhood of d J in C in which R{x) = 0(e £S ). 
Proof. The function 

7](x) := I (R(x) - e~ isx - q(x)) e isG{x) \ 2 (33) 

is bounded and subharmonic in Q and it follows from Lemmas 6 and 7 that its 
boundary function satisfies 

rj(x) \dx\ = 0(s^ 1 ). 



Hence, from general considerations, rj(x) = 0(s x ) uniformly on compact subsets of Q. 
Combining this with the estimate of Lemma 8(a) proves our first assertion. (Observe 
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that the error term in the statement is exponentially large and so dominates the term 
e~ isx .) To prove the second we replace G(x) in (|33"D by the analogous function where 
J is replace by J minus small neighborhoods of its boundary points. We then get an 
estimate for R(x) — e~ tsx — q(x) also in a neighborhood of d J, apply Lemma 9(b), 
and deduce the second assertion. □ 



IV. Demonstration of (|j) 



Let us first evaluate asymptotically the integral 

X:= f \h+(x) +h-(x) -e~ lsx \ 2 dx. (34) 

JE 

It can be written as 



E 



h + (x)(h+(x) -e' lsx )dx + / h-(x)(h-(x) -e~ lsx )dx 



(h+(x) e lsx + hJx) e lsx -l)dx+ / (h + (x) h_{x) + h_(x) h+(x)) dx. 

JE 

Since this is real it equals the real part of what is obtained by replacing the third 
integral by its complex conjugate, so the above equals 



Re |j& h(x) (h(x) - e~ tsx ) dx - J E (h+(x) e~ lsx + h_(x) e~ tsx - 1) dx 
+2 Re / h+(x) h-(x) dx = Re (2i + J 2 ), 

JE 

say. Thus I = Re (Ii +X 2 ). We shall evaluate T\ exactly and show that X 2 is 0(s _1 ). 
Evaluation ofX\. 

We write our integrals over E as limits as r — > oo of the corresponding integrals 
over E fl [— r, r], which can be combined. We obtain 

2i = r lhn | j h(x)(h(x) - 2 e~ isx ) dx + 2r| - | J|. (35) 

Now we can almost, but not quite, apply Lemma 4 here with 

g{x) = h{x)-2e~ isx . 

The problem is that this function does not satisfy the first hypothesis of the lemma, 
but we shall just modify its proof. Recall the notation there; for convenience we drop 
all subscripts T. We have 

h(x) (h(x) - 2 e~ lsx ) dx=l j(x) e isG{x) (e~ isG{x) k(x) - 2 e~ isx ) dx 

L JE 
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_j(x)(k(x) -2e ls{G{x) - x) )dx, 

where C + and C~ are semi-circles described clockwise in the upper and lower half- 
planes, respectively, joining — r and r. Since j(x) and k(x) are 0(x _1 ) as x j oo and 
Im (G(x) — x) > there, the integral over C~ tends to zero as r — > oo. As for C + , 
both j(x) and tend to 1 as x f oo, and are even analytic there. Suppose 

j(x) = 1 + - + 0(1), fe(z) = 1 + — + 0(\) as x too. 



If we use (|15|) then we find that the integrand equals 

-1 + ^ + ^^ + 0(1) 

XXX 2 

as x f oo. Hence the integral over C + equals — 2r + 7r (2a 2 s + i(j — fc )) + o(l) as 
r — > oo and (^5|) gives 

Ji = 2na 2 s + 7ri ( jo - k ) — \ J\, 

whence 

Re Ii = 2na 2 s + tt Im (A; — jo) — 1^1- 

Observe that 

fc = limx(k(x) — 1) = limx \ogk(x), 
and similarly for j(x). Thus 

k(x) 



Im (fco — jo) = hm a; arg- 



Now j(x) is obtained from k(x) by changing the signs of the in ([H]), so that 

•7W fe=o 

Recalling the normalizations fl2"2"|) we deduce that 



m—l 

Im (/c - jo) = Y] £k lim x#(x, x k ), 

• — ' x— >oo+0j 
fe=0 

and so 

m— 1 

Re Xi + | J| = 27ra2S + 7T V" £fc lim xg(x, Xk). (36) 

k=0 



Proof that X 2 = O^ 1 ): 
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Since Gq(x) = —Go(x) and G(x) = Gq{x) — a\ we have, with an obvious notation, 

=k+{x)k-{x) e 2i < G +^ +a ^. 

Again we use the facts that the product of the exponential Green and Neumann 
functions in (pi] ) is smooth, uniformly in s, and that the local behavior of the other 
functions arising here is determined locally. In particular k + (x) k-(x) is 0(|x| -1 ) at 
infinity, in fact for each n its n'th derivative is 0(|x|~ n ), and it has the behavior of 
(x — afc)~2 near and analogous behavior near fa. The function G+(x) has the 
behavior of x near infinity by fll5|), and 

G'(x) = P{X) 



q{x) 

which is positive on E + . It follows easily from these facts that if <p(x) is any C a 
function with compact support which is identically 1 on a neghborhood of dJ then 

J {I- tp{x))k^{x)k-{x) e 2isG+ix) dx = 0(s~ n ) 

for any n, while 

' (p(x)k^x) k-{x) e 2lsG+{x) dx 



ip(x) — =r- dx, 
o Jx 



has the behavior of 



where ip is a C°° function with compact support. Such an integral is 0(s _1 ). 
Recapitulation: 

We have shown that X, the integral in ([34]), equals 

m—l 



27ia 2 s + 7i Ek lim xg(x, Xk) — \ J\ + 0(s ] 

f— ' x— >oo+0i 
k=0 

In particular X = O(s) and so from Lemma 7 we get 

/ \R(x) -e isx \ 2 dx = J + 0(1), 
Je 

and this and (0) give 

-7r^-log det (I-K s ) =X+ \J\+0(1). 

ds 
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If the error term 0(1) here were 0(s *) we would have established (^) with 

m— 1 

c 1 = 2a 2 , c 2 (s) = V £*. lim xg(x, x k ). (37) 

f— * x— ►oo+Oj 
«=0 

As it stands, of course, we have only proved the weaker statement 

-4- log det (J - K.) = as + 0(1). 

as 

Conclusion: 

Here is how to get the stronger statement. The details of the analogous argument 
for orthogonal polynomials can be found in ]T5| . 
It follows from ([]) that 

| J g{x) e isx dx\ <\\g\\ \\ R(x) - e~ isx \\ 

for all g G S s , the norms being that of the space L 2 (E), and that this becomes an 
equality when g(x) = R(x) — e~ tsx . From this and (|7[) we deduce 

-4 log det (I - K.) = \J\ + max L&g^£^l! , 
as ' ges s \\ g \\ z 

Denote the maximum on the right by Mq. It follows from Lemma 4, in the form 

/ (h+(x) + hjx) - e~ isx ) g(x) alx = f g(x) e isx dx 
Je Jj 

(first for Fourier transforms of smooth functions in L 2 {— s, s) and then by continuity 
for all such functions), that 

M < [ \h + (x) + h^{x)-e- isx \ 2 dx. 

JE 

The right side is just what we have called X and so we know that 

M < 7T Cl + ttc 2 (s) - \J\ +0(s- 1 ). (38) 

Now we modify the the extremal problem by replacing || g \\ by the norm of g in 
L 2 (E, ws) where the weight function wg(x) is equal to 1 except on the 2£-neighborhood 
of dE, equal to q(x)~2 in the (^-neighborhood of dE and, say, linear in between. 
Denote by the corresponding extremum for the weight function w$. Clearly 

M S <M . (39) 
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There exist reproducing functions analogous to k-p in ( pT]) for any weight function 
satisfying quite general conditions, so that \dx\ in the statement of Lemma 3 can be 
replaced by ws{x) \dx\. The extra ingredient is the introduction of a certain nonzero 
function in Ti whose absolute value on E equals w$- We then proceed exactly as before, 
defining the analogous h s by (p3|) and q by (|27|). The main point is that because ws 
has the behavior of q(x)~z near dE the corresponding reproducing function kr has 
the behavior of q(x)z, with the result that the factor t 1//4 does not appear in the 
denominator in (|32j). This improves the error estimate in Lemma 6 to 0(s~ 2 ) in this 
case, and so in the end we obtain 



M s = vrci + vrc 2 , 5 (s) - \J\ + Ois- 1 ] 



where 02,5(5) is the function in ([37]) associated with the weight function ws- Because 
of the continuity of this function in the class T as well as the weight function w, we 
deduce that for any e > there is a 5 such that 

M 5 > vrci + vrc 2 (s) -\J\- e + o(l). 

Putting this together with fl38|) and (|39| ) gives the desired result. 



When J is a single interval of length 1, say J = [— |, |], then G(a; 



2 , 2J , ^ w ^ 4 , 



and we find that ci = 1/4, c 2 = 0. Thus (§) in this case is (0) with the improved 
error term 0(s _1 ). Notice that this is consistent with ([I]) since the first power of s 
does not appear in the exponent. 

It is easy to see that c 2 (s) is periodic when m = 2 and to compute its period. For 
its value depends only on the class of e isG ^ x \ which in turn depends only on the value 
modulo 2tc of s times the quantity A^q given by fll3|). Hence the period of c 2 (s) 
equals 

1 f a ? p(x) 



7T 



Pi sjqjx] 



dx } . (40) 



V. Demonstration of (|5D 

It follows from Lemma 9(b) that for each e > the contribution to the integral in 
(|8|) of a sufficiently small neighborhood of d J is 0(e es ). The main contribution will 
be from certain interior points of J, as we shall see. Outside any neighborhood of dJ 
the asymptotics of R(x) are given by Lemma 9(a), which may be rewritten as 

R(x) = e~ isG{x) (ife(x) + 0(s~2)). 
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It follows from general considerations that the derivative of the O term above is also 
Ofs-5). We use the formula 



u' 



(uv)' uv = \u\ ( — \v\ + v'v) 



u 



with u equal to the exponential factor above and v the other to write 

R'(x) R{x) = e 2slmG{x) {-tsG'(x) \k(x) + 0{s^)\ 2 + k'(x) k(x) + O(s^)} . 

Since G'(x) is purely imaginary in J the first term in the braces is purely real, and so 
does not contribute to the right side of @, and we are left with the computation of 



3 2slm G(x) 



(Im k'(x) k(x) + 0(s-*))dx. (41) 



The integration is taken over J with its little neighborhood of dJ removed. The 
main contribution to the integral will come from the point or points where Im G(x) 
achieves its maximum, and these will be among the zeros of p(x). Denote these zeros 
by Zi, ■ ■ ■ , z m , one lying in each interval of J. Since at each Z{ 

d 2 t r>t \ T d P( X ) \p'( X )\ , An \ 

— Im G(x) = Im ^—7= = - r—- , 42 



we find by standard asymptotics that (f41|) equals 

flf e2 ,i m a M {M^Im k'{ Zi ) k[z7) + O(s^) 
V s I \p'{Zi)\2 



i=l 

We have from ( PT|) 



■y m— 1 



= 2 II l*(*i> X k *X k )$(Zi, x k ) £k \ 



k=l 



and, using in addition to fl2~l~|) the very last part of (|4"2|), we find that 



Combining these last two formulas gives Im k'(zi) k(zi). Comparing with (|8]) shows 
that (|5|) is established, with 

c 3 = 2 maxfc Im G(z k ) 
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and 



C4(s) -40Ftib^)|^ 



\p'(z-)\ m ~ l ~ 

' "' J2( N ( z i> x k*Xk) +e k g{z h x k *x k )) 



\q(zi)\ 



k=l 



x [] x k *x k )^(z h x k ) £k \\, (43) 

k=l ) 

where the outer sum is taken over those % for which 

Im G(zi) = maxfc Im G(z k ). 

There is an explicit representation of the quantities Im G(z k ). It follows from flTT|) 
and ( |T2"D that for z G (a k , (3 k ) 

ImG(z) = f Z ^EkLdx, 



° k y/\q(x)\ 



the sign being that of p(a k ). The maximum of this, its value at z k , is the integral to 
z k , while the integral from z k to (3 k is the negative of this since the sum of the two is 
zero. Hence 

Im G(, fe ) = i T J^MUz, 



2 

and so 



a I, 



q(x)\ 



f Pk \p(x)\ , 
C3 = maxfc / ax. 



\q(x)\ 



When J is a single interval of length 1 this equals 1. 
When m = 1 (^) becomes 

1 \P'(^ 
4^\q(z)\y^ 



simply, and when the length of J is 1 this is l/2v2vrs and we recover the formula 
When m = 2 the function c 4 (s) is periodic with the same period (|^) as c 2 (s). 
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